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D ■ Abstract. Duality relations between continuous-state and discrete-state stochastic 

processes with continuous-time have already been studied and used in various research 
fields. We propose extended duality relations, which enable us to derive discrete-state 
stochastic processes from arbitrary diffusion-type partial differential equations. The 

& \ derivation is based on the Doi-Peliti formalism and the algebraic probability theory, 

"£i ■ and it will be clarified that additional states for the discrete-state stochastic processes 

r^ , must be considered in some cases. 
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1. Introduction 

Duality is one of the widely used techniques to investigate interacting particle systems 
(for example, see [1]). Using the duality concept, a stochastic process may be 
connected to a different kind of stochastic processes; for example, duality relations 
between a continuous-state stochastic process (a stochastic differential equation) and 
a discrete-state stochastic process (a birth-death process) have been investigated a lot 
in various research contexts ranging from population genetics [2|,|3] to nonequilibrium 
heat-conduction problems [I]. 

Although the constructions of a dual process have mainly been performed with 
a kind of heuristic ways, recently general schemes to derive the dual process have 
been investigated [HIE]. In [6], the Doi-Peliti formalism [7H9] has been used to derive 
the duality relations between continuous-state and discrete-state stochastic processes. 
As a consequence, it is possible to derive partial differential equations from arbitrary 
birth-death processes describing chemical reactions or ecological systems. Note that 
for a second-order partial differential equation it is sometimes possible to derive a 
corresponding stochastic differential equation, and hence we can derive the duality 
relations between the continuous-state and discrete-state stochastic processes. However, 
the reverse direction has been unclear. That is, it was unclear whether birth-death 
processes can be derived from from 'arbitrary' partial differential equations or not; 
the Doi-Peliti formalism can derive a kind of linear operators from partial differential 
equations, but there is no guarantee that the derived operator can be interpreted as 
that of birth-death processes. 

In the present paper, we propose extensions of the duality concept, which enables 
us to derive a birth-death process from an arbitrary diffusion-type partial differential 
equation. In order to connect the partial differential equation with a discrete-state 
'stochastic' process, additional states should be introduced in some cases. In order 
to discuss such additional states, we employ the scheme of the algebraic probability 
theory. Using the extended duality relations, it is possible to evaluate moments of the 
continuous-state stochastic processes or partial differential equations from the solutions 
of the discrete-state stochastic processes. 

The construction of the present paper is as follows. In section 2, the conventional 
duality is briefly reviewed, in which the derivation of the dual process via the Doi- 
Peliti formalism is shortly explained. Section 3 gives an extension of the duality, i.e., 
a weight factor is introduced. The main claim of the present paper is discussed in 
section 4, in which not only a weight factor, but also additional states are introduced 
in order to define the extended duality relations. Section 5 gives discussions. Although 
all discussions in the present paper will be performed by using some specific examples 
for clarity, it is straightforward to apply the discussions for general cases. 
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2. Conventional duality 

In this section, we briefly explain the conventional duality relations between continuous- 
state and discrete-state stochastic processes. The discussions and derivations are 
basically based on [6]. In [6], it has been pointed out that the correspondence between 
the generating function approach and the Doi-Peliti formalism is available to derive the 
duality. However, we noticed that there is no need to use the correspondence explicitly. 
We here give the derivation without using the generating function approach. 

2.1. Definition of the duality 

The conventional duality between a continuous-state stochastic process and a birth- 
death process is defined as follows. Suppose that (zt)t>o G R is a continuous-state 
and continuous-time Markov process, and that (n t )t>o G N is a birth-death process. 
Let E n denote the expectation with respect to the process (n t ) t > which starts from 
n . The process (n t )t>o is said to be dual to (z t )t>o with respect to a duality function 
D : R x N -> R if for all z E R, n G N and t > we have 

E n [D{z ,n t )]=E z [D{z t ,n )], (1) 

where E^ is expectation in the process (z t )t>o starting from z . 

2.2. Problem settings 

In order to see the duality relations in more explicit forms, we here consider a concrete 
partial differential equation, which is related to the stochastic Fisher and Kolmogorov- 
Petrovsky-Piscounov (sFKPP) equation; the sFKPP equation plays an important role 
in the study of the front-propagation problems [T0HT3] . and even in a QCD context |14j . 
The partial differential equation with constant parameters 7 and a, 

§£P(z, t) = -j^ [-7*(1 - z)p{z, t)\ + \^ y^H - z)p(z, t)] , (2) 

can be interpreted as the following stochastic differential equation: 



dz = - 7 z(l - z)dt + ay/z(l - z) 6W(t), (3) 

where W(t) is a Wiener process. It has been shown that the n-th moment of the 
stochastic differential equation in ([3]) or the partial differential equation in (J5J), defined 
as 

/<x> 
dzp(z,t)z n , (4) 

■00 

can be evaluated without solving these equations directly. That is, the evaluation is 

performed from a corresponding birth-death process via the duality relations introduced 

in (HJ. 
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2. 3. Derivation of the dual process and duality function 

The problem here is: how can we find the corresponding dual process and the duality 
function? In order to find them in a unified way, we first introduce creation and 
annihilation operators in order to discuss stochastic processes [7H9]. This method is 
called the Doi-Peliti formalism, the second quantization method, or the field-theoretic 
method, and so on, and it has been used to investigate stochastic processes mainly 
from the view point of perturbation calculations or renormalization group methods (for 
example, see [16]). 

In the Doi-Peliti formalism, the following bosonic commutation relation is 
introduced: 

[a, a'] = <W — a^a = 1, [a, a] = [a', a*] = 0, (5) 

where a^ is the creation operator, and a the annihilation operator. Each operator works 
on a ket vector in Fock space \n) as follows: 

a'\n) = \n + 1), a\n) = n\n — 1), (6) 

and the vacuum state |0) is characterized by o|0) = 0. The bra vector (m\ is naturally 
introduced from the inner product defined as 

(rn\n) = <5 m ,„,n!, (7) 

where 5 m , n is the Kronecker delta. 

In order to connect the partial differential equation in (J2J) with the discussions 
in terms of the Doi-Peliti formalism, it is convenient to introduce the following linear 
operators: 

d \ da 2 d 2 

;. S J = -7;(l--% + T -~(l-%j (8) 

and 

i.e., L* is the adjoint operator of L. Hence, the partial differential equation in ((2]) can 

be rewritten as follows: 

d ( d\ 

-p(z,t) = L*U — )p(z,t). (10) 

Next, the following state vector (p(t)\ characterized by the probability density p(z,t) is 
introduced: 



(p(t)\= dzp(z,t)(z\, (11) 

J — oo 

where (z\ is a coherent state of a\ i.e., 

(z\ = <0|e- (12) 

which satisfies 

(z\a) = z{z\, (13) 
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and z is assumed to be a real variable. 

We here construct the linear operator L {a), a) by simply replacing z and J^ with 
a' and a, respectively; 

2 

L( a t, a) = -7a f (l - a f )a + —a) {I - o f )o 2 . (14) 

Note that the linear operator L(a^,a) is written in the normal order; all annihilation 
operators a are to the right of all creation operators in the products. Then, using the 
following formula, 

{z\{a)) k a l = (z\z k a l = z k (0\ (J-) e za = z k (J-\ (z\, (15) 

we have the following identity: 



(p(t)\L [a\ a) = / dzp(z, t)(z\L (a^, a) 

J — oo 

= / dzp(z,t)L (z, — J {z\L{a) ) a) 

(z\. (16) 



dz 



oo 



Z, dz~} 



Hence, the solution of the partial differential equation in ([2]) is obtained formally from 
the following state vector: 

(p(t)| = <p(0)|exp[L(ot,o)t]. (17) 

Note that it is also possible to consider cases with time-dependent coefficients. For 
simplicity, we here restrict ourselves to the case with time-independent coefficients. 

Next, a state vector spanned by the basis of the Doi-Peliti formalism is introduced 
as follows: 



\P(t)) = ^P(n,t)\n), (18) 

n=0 

where we assume that P(n,t) is adequately normalised as 

oo 

£p(M) = i. (19) 

n=0 

Assuming that the action of the linear operator L(a\ a) to the state vector \P(t)) gives 
the time-evolution of the state vector, i.e., 

exp [L(a\a)t] |P(0)) = \P(t)), (20) 

the following duality is obtained: 

(P(0)|P(*)> = {p(0)|exp[L(at,o)t]|P(0)) = (p(t)|P(0)). (21) 

Hence, the duality in (I21p is explicitly written as 

/oo °° /*oo °° 

dz^p(z,0)P{n,t)z n = / dzY,v{z,t)P{n,d)z n , (22) 

00 n=0 ^~°° n=0 



Extended duality relations between birth-death processes and PDE 6 

where we used the following identity: 

oo 1 oo . oo 1 

W = (°l E ^"^^ = E ^ m H^> = E -i* m Sm,nm\ = z\ (23) 



mi * — ' mi * — ' mi 

m=0 m=0 m=0 



And, therefore, we conclude that the duality function is given as 

D{z,n) = z n . (24) 

If we set the initial conditions for n and z as a Kronecker delta function and a Dirac 
delta function respectively, the duality relation in (Q]) is recovered; using the explicit 
expression for the duality function, we have the following duality relation: 

E n [z^] = E z [z?°] . (25) 

The remaining problem in the above discussions is the verification of (|20|) . Is it 
possible to interpret ( |20l) as the time-evolution of the probability distribution P(n,t)l 
If not, f l20|) is not valid and the duality relation is not derived. 

In order to check the validity, we split the time-evolution operator exp(L(a>, a)t) 
as a product of exp(L(a^, a) At) with small At. The expansion of exp(L(a^, a) At) gives 

exp(L(a f , a) At) ~ 1 + L(a\ a) At 



2 

1 -7a f aAt- — (a)) 2 a 2 At 



2 

7 (a f ) 2 aAt + — a j a 2 At, (26) 



where 1 is the identity operator. Noting a^a\n) = n\n), the expansion in ( 1261) can be 
interpreted as follows. Assume that the current state is |n). The second term 7(a^) 2 aAt 
means that there is a state change n — > n + 1 with probability ^nAt. The third term 
(o- 2 /2)a^a 2 At corresponds to a change n — > n — 1 with probability (a 2 /2)n(n — l)At. 
The first term corresponds to the case in which no transition occurs. The sum of the 
probabilities for these events is equal to one, and hence the operator exp(L(a^, a) At) 
gives the time-evolution for the probability distribution. More explicitly, we consider 
the following time evolution equation for the state vector \P(t)): 

j t \P{t)) = L{a\a)\P{t)), (27) 

which gives the following master equation for P(n,t) by comparing the coefficients of 

\n) of the left and right hand sides in (127)1 : 

d 

—P(n,t) = 7(n - l)P(n - l,t) - jnP{n,t) 

(n + l)n , . n(n — 1) , , . 

+ o- 2 V 2 ; P(n + 1, t) - a 2 K 2 } P(n, t). (28) 

The master equation in (|28|) suggests the following birth-coagulation process for particles 
A: 

Reaction 1: A — > A + A, 

Reaction 2: A + A -> A, 



i.e., 



n — > n + 1 at rate '-/n, 

n — > n — 1 at rate a 2 n{n — l)/2, 
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where n is the number of particles A. Hence, we can verify (j20p . and the partial 
differential equation in (J2J) or the continuous stochastic process in ([3]) is connected to 
the birth-death process in ( f28l) via the duality relation in ( |25l) . 

In general, the discrete-state stochastic processes are easily solved compared with 
the continuous-state stochastic processes, and actually the duality relations have been 
mainly used to evaluate various quantities for the continuous-state stochastic processes 
(e.g., see 0). 

3. Duality with weight factors 

3.1. Problem settings 

Instead of the original problem in (J2J), we next consider the following partial differential 
equation: 

dt P ^ f) = ~~d~z [_7Z(1 ~ Z)P ^ t)] + 2d^ ^^^ ^ ■ (29) 

Note that the coefficient in the second term in the right hand side is different from 
(T5]). The partial differential equation can be interpreted as the following stochastic 
differential equation: 

dz = -jz(l-z)dt + a^dW(t). (30) 

The corresponding adjoint operator L* (z, Jj) is 

_j = __ [ - 7z(1 - 2)1 + __^, (31) 

and the linear operator in terms of the creation and annihilation operators is 

L{a\a) = -7a t (l-a t )a + — a) a 2 . (32) 

One may expect that the linear operator L(a\a) describes time-evolutions of a 
certain type of stochastic process. However, the same discussions in section 12.31 do not 
give a stochastic process; the probability conservation is violated and (1201) cannot be 
verified. In order to derive a dual stochastic process, it is needed to extend the definition 
of the duality using a weight factor. 

We here note that the weight factor has already been introduced for the duality 
relations in spin systems [I]; we apply it to the duality between continuous-state and 
discrete-state stochastic processes via the Doi-Peliti method. 

3.2. Duality 

The second term in ( |32|) changes the states; it contains a^a 2 , and then there is a state 
change n — > n — 1. However, there is no corresponding term for the probability 
conservation. Hence, a simple solution for this problem is to add an additional 



Extended duality relations between birth-death processes and PDE 8 

term, which does not change the state and complements the lack of the probability 
conservation. That is, we rewrite the linear operator L{a\a) in (1321) as 

2 2 

L{a\ a) = -7^(1 - a))a + — ^(1 - a^a 2 + — (a f ) V 

= L'(a\a) + V{a)a), (33) 

9 

L'(a\ a) = -7a f (l - a ] )a + — a f (l - a))a 2 (34) 



where 



and 

2 

y( fl t a ) = ^_ [( a t a )2 _ a t a ] . (35) 

Note that the operator ^(0^0) is expressed only in terms of a)a, and therefore the 
operator V(a^a) does not change the state \n). Considering a small time interval At, 
we have 

exp [L(a\ a) At] \n) = exp [V(a^a)At + L'(a\ a) At] \n) 
~ (1 + y(a f a)At + L'(a), a) At) \n) 
= {l + V(n)At + L\a),a)At)\n), (36) 

and we can replace exp[V(a^a)At + L'(a\ a) At] with exp[V(n)At + L'(a\a)At] for 
the small time interval At, which makes exp(V(n)At) a scalar value. Hence, only 
the factor exp[L'(cr , a) At] plays a role as a time-evolution operator. The factor 
exp \U (a', a) At] can be interpreted as the time-evolution operator for the birth-death 
process in section 12.31 and we have 

exp [L\a),a)t] |P(0)) = \P{t)), (37) 

instead of ( 120]) . As a consequence, the duality relation should be modified a little, and 
finally the following duality relation is obtained instead of ( 125]) : 



E„ 



exp J J V{n t ,)dt'Xz^ 



E z [z?>] . (31 



Note that the expectation on the left hand side of (1381) should be taken for all the paths 



because the weight factor exp(/„ V(n t >)dt') depends on the path. 



4. Duality with weight factors and additional states 

4-1. Problem settings 

In this section, we deal with more complicated cases, which is seen in the following 
example: 

J^(*> *) = ~§- z [7*(1 - *)p(z, t)l + l§^ OM* " *M*. *)] » ( 39 ) 

which can be interpreted as 

dz = 72(1 - z)dt + ay/z(l - z) dW(t). (40) 
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The adjoint operator L* (z, -SA is written as 

''(^I'-sW-^^-'^ (41) 

and in terms of the creation and annihilation operators, L(a\ a) becomes 

2 

L(a\ a) = 7a f (l - a))a + ^-a f (l - a f )a 2 

2 

= 7a t a-7(a t ) 2 a + — a f (l -a f )a 2 . (42) 

Notice that there is the term, — ^(a^) 2 a, in (142]) . This term changes the state as 
n — > n+1, but the transition probability is negative. (In order to see it, use the expansion 
similar to (T26]) .) It is impossible to consider such negative transition probabilities for 
usual stochastic processes. In addition, we cannot use the similar techniques introduced 
in section 13. 2( there is no way to cancel the effects of the transition with the negative 
probability and to split the operators in a simple form as (136]) . 

In the following discussions, we introduce additional operators and an additional 
state, based on the algebraic probability theory. This additional state enables us to deal 
with the transitions with the negative probability adequately. 

Jf.,2. Additional states and duality 

We here give a brief review of the algebraic probability. The algebraic probability is also 
called quantum probability or noncommutative probability, and, for example, see [17] 
for details. 

The algebraic probability is based on a *-algebra and a state defined on the *- 
algebra. In order to define a *-algebra, an involution is introduced as follows: the 
involution is defined as a map a H- a* on an algebra A and satisfies 

(a + b)* = a* + b\ (Aa)* = Aa*, (ab)* = b*a\ (a*)* = a (43) 

for a, b G A and A 6 C. A *-algebra is an algebra equipped with an involution. Next, we 
define a linear function ip on a *-algebra A with values in C, which satisfies the following 
two properties: (i) (positivity) tp(a*a) > 0; (ii) (normality) <p(1a) = 1, where 1^ is the 
identity. A linear function ip with the positivity and normality is called a state, and an 
algebraic probability space is a pair (A, tp) of a *-algebra A and a state tp on it. We here 
note that the Doi-Peliti formalism is also one of examples of the algebraic probability 
space [T8] . 

In order to discuss the negative transition probability, it is useful to define the 
following algebraic probability space. New state vectors, |+), |— ), (+|, and (— |, are 
introduced, which satisfy 

(+1+) = (-|-) = i, <+!-> = H+> = 0. (44) 

Next, new operators b = b* and iV = N* are defined as 

6|+> = H, 6H = l+>, (+\b=(-l (-\b=(+\, (45) 

iV|+) = |+}, TV]— > = — 1— >, (+\N=(N\, (-\N =(-|-. (46) 
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Note that the operator N plays a role as a number operator. We consider an *-algebra 
A constructed by 0, 1, b = b*, and N = N*. Using a state vector \p(t)) defined as 

\p(t))=p(+,t)\+)+p(-t)\-), (47) 

where p(+, t) and p(— , t) are probabilities with which the state is in + or — , respectively, 
a state 4>(a>) for a G A is defined as 

0(«) = ((+! + (-!) a|p(t)>, (48) 

for any a£ A Noting that p(+,t) +p(—,t) = 1, it is straightforward to confirm that 
the above construction satisfies the conditions of the algebraic probability space. 

Using the above additional states |+) and |— ), we consider the following combined 
stochastic process. 

Firstly, a new state vector \ip(t)) for the birth-death process is introduced as 

oo 

M*)> = £ E P(n,s,t)\n,s), (49) 

n=0 s={+ ,-} 

where P(n, s, t) is the probability with which the number of particles in the birth-death 
process is n and the sign of the state is s G {+, — } at time t, and \n, s) = \n) <8> \s). 
Note that the operators a' and a act only on the particle numbers n, and the operator 
b affects only on the sign of the state s G {+,—}• That is, 

a'|n, s) — \n+ 1, s), a\n,s)—n\n — l,s), (50) 

for s G {+, — }, and 

b\n,+) = \n,-), b\n,-) = \n,+). (51) 

In addition, we define (p(t)\ as follows: 

/oo 
dzp(z,t)[(z,+\-(z,-\]. (52) 

-oo 

Using the above definition, there is no need to define the sign of the states, + or — , for 
the partial differential equations, as we will see soon. 

Secondly, we construct a linear operator acting on the new state vectors \ip(t)) 
based on the linear operator L(a\ a) in (142]) . as follows: 

2 

L(a\ a, b) = ^a^a + b^(a)) 2 a + —a) (I — a))a 2 

2 

= 27a f a - 7a f (l - ba))a -\ a){\ - a ] )a 2 

= L'(a\a,b) + V(a^a), (53) 



where 



and 



2 

L'(a\ a, b) = - 7 a f (l - ba))a + — a f (l - a r )a 2 (54) 



V(a^a) = 2 7 a f a. (55) 
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Note that for the negative sign in the second term in (j42p is replaced with the new 
operator b. 

Thirdly, we check the action of L{a\ a, b) for (p(t)\; 

(p{t)\L{a\a,b) 

dzp(z, t) [{z, +| - (z,—\] L(a\ a, b) 



f d d d 

dZ I Q~\~ 1ZP ( Z ^ ^' + ' ~~Q~Z \- 1Z P ( Z,t ^ ^'~l + Q-b z P( Z ' t )] (*'~l 

+ ^ [ 7 ^p(«, t)] (z, -| + 1^ [a 2 (l - z)zp(z, t)] [(z, +| - (z, -|] } . (56) 
On the other hands 



j - 



A(p( f )| = jT d8 (^P(-M)) [<*,+! ~ (Z,-W • 



(57) 



Comparing the coefficients of (z,+\ (or (z, — |) in (156]) and f )57|) . we have the original 
partial differential equation in ( 1391) . Hence, the new time-evolution operator defined in 
fl53|) adequately recovers the original problem. 

Fourthly, we investigate the action of the time-evolution operator L\a* , a, 6) in ( 1531) . 
(The effects of V(a^a) are the same as in section |372| i.e., a weight factor.) After some 
calculations, the following master equations for P(n, s, t) are obtained: 

d 
P(n, +, t) = 7(n - l)P(n - 1, -, t) - -ynP(n, +, t) 



dt 



(n+l)n„, „ . 9 n(n — 1) „, , ,_ 

+ o- 2 - — -^P(n + l,+,t)-a 2 ^- — '-P(n,+,t), (51 



—P(n, -, i) = 7 (n - l)P(n - 1, +, t) - 7 n?(n, -, t) 

+ ^ 2 ^ P(n + 1, -, t) - a 2 y 2 } P(n, -, t). (59) 

The master equations in (|58|) and (159]) can be interpreted as interacting particle systems 
with internal states, i.e., the sign of the states; the birth-death process described by (158]) 
and (1591) becomes as follows: 



Reaction 1: A — > A + A and change the sign of the states, 
Reaction 2: A + A -> A, 



n — > n + 1 and "+ state — > — state" (or "— state — > + state") at rate 7 n, 

n — ?• n — 1, at rate o~ 2 n(n — l)/2, 

where n is the number of particles A. Only when the reaction 1 occurs, the sign of the 
state (the internal state) is changed. 

Finally, the duality relation for the above extended stochastic processes is written 
as follows: 

f /"' 



expj f VMdAz? -E n ,_ expj f V{n v )dt\z^ 



E z [z?°] , (60) 
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where E n+ and E„ 5 _ mean the expectations with respect to P(n,+,t) and P(n,—,t), 
respectively. In addition, as an initial state, we here set p(z, 0) = 5(z — z ) for the partial 
differential equation and P(tiq, +, 0) = 1 for the discrete-state stochastic process. We 
note that it is possible to extend the duality relations for arbitrary initial conditions. 

5. Discussions 

In the present paper, the duality concepts were extended for various cases by using 
some specific examples. It is straightforward to apply the scheme to multivariate cases. 
And when we consider the multivariate cases, the additional state in section 4 can be 
reinterpreted in another way. That is, we introduce a new variable y, which obeys the 
time-evolution with 

4- a («) 

When we set the initial state of y as y(0) = y$ = — 1, f|40l) can be rewritten as follows: 

dz = (~/z + -?yz 2 ) dt + ay/z(l - z)dW{t), (62) 

because y is always — 1. In the dual process, y is replaced with the creation operator 

at, which creates an additional particle corresponding to y. Denoting the dual discrete 

variables for z and y at time t as n\ and rtf respectively, the duality relation is written 

as 

"t -\ i ,- 

(63) 



expj^ V(n z t ,)dt^zfyf 



^zo,3/0 



ZtVt 



{.Jo J 

Since y t = yo = — 1, we can recover ( 1601) by setting n\ = 0. 

The scheme is also applicable for equations with more complicated functions. For 
example, we consider the following stochastic differential equation: 

ds = sin(s)dt + dVF(t). (64) 

In order to construct the corresponding dual process, one might consider that the 
function sin(s) should be replaced with sin(aj). However, it is not easy to interpret 
this term within the dual 'stochastic' process. Instead, we introduce two additional 
variables: 

u = sin(s), v = cos(s). (65) 

(66) 
(67) 

(68) 

and the variables s, u, and v are interpreted as a\, a* , and aj in the corresponding dual 
process. That is, by using the discrete stochastic process with three variables, nf, n", 
n v t , we can construct the dual process for (J51]) . 



Hence, we have 




ds 
dt 


dW(t) 

= '"+ * ■ 


du 


ds ds 


dt 


= cos(s)-— = v—, 
K ' dt dt 


dv 


ds ds 


dt 


= -sin s — = — W— 
v J dt dt 
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It may be difficult to solve derived discrete-state stochastic processes analytically. 
Even in the cases, the extended duality concepts are useful. For example, when one 
wants to obtain some moments for a stochastic differential equation numerically, a kind 
of approximation scheme is needed. The basic one is the Euler-Maruyama scheme [19] , in 
which a time-discretization is inevitable. Hence, in order to make good approximations, 
it is necessary to choose an adequate time-discretization. In contrast, there is no need 
to use such time-discretization for simulations of the derived discrete-state stochastic 
process; there are various simulation schemes for birth-death processes without using 
the time-discretization. One of the famous algorithms is the Gillespie algorithm [20] , in 
which time-intervals between events are selected directly from exponential distributions, 
so that the time-discretization is not needed. In addition, the Gillespie algorithm is not 
an approximate one, hence rapid and precise calculations could be available by using 
the dual stochastic process. 
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